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Abstract 

In this paper we investigate the Bose-Einstein condensation of massive spin-1 
particles in an Einstein universe. The system is considered under relativistic 
conditions taking into consideration the possibility of particle-antiparticle pair 
production. An exact expression for the charge density is obtained, then certain 
approximations are employed in order to obtain the solutions in closed form. 
A discussion of the approximations employed in this and other work is given. 
The effects of finite-size and spin-curvature coupling are emphasized. 
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1 Introduction 



Since the early work of Altaie (1978), the study of Bose-Einstein condensation (BEC) 
in curved space has attracted the interest of many authors (Carvalho and Rosa 1980, 
Singh and Pathria 1984, Parker and Zhang 1991, Smith and Toms 1996, and Trucks 
1998). This interest comes in the context of trying to understand the thermodynamics 
of the early universe and the role played by the finiteness of the space in determining 
the thermal behavior of bosons in the respective systems. 

Most studies employ the static Einstein universe as the underlying geometry for 
the system, since the thermodynamic equilibrium in this universe can be defined 
without ambiguity (Altaie and Dowker 1978). In a time developing spacetime no such 
luxury is enjoyed, except for the Robertson- Walker spacetime which is conformally 
static (Kennedy 1978). 

In an earlier work (Altaie 1978), we considered the cases of non-relativistic BEC 
of spin and spin 1 bosons in an Einstein universe. We found that the finiteness 
of the system resulted in "smoothing out" the singularities of the thermodynamic 
functions found in the infinite systems, the enhancement of the condensate fraction 
and the displacement of the specific-heat maximum toward higher temperatures. 

Carvalho and Rosa (1980) considered the case of relativistic scalar field in an 
Einstein universe in an effort to show that the finite size effects are negligible in 
comparison with the relativistic effects. However, a close look at their work shows 
that this comparison was done at implicitly large value of the radius a of the Einstein 
universe. And since finite size effects are proportional to 1/a the finite size effects were 
found to be negligible. Moreover, the consideration of relativistic and ultra-relativistic 
BEC should take into account the possibility of particle-antiparticle production, since 
at temperatures greater than the rest mass of the particles, quantum field theory 
requires the inclusion of such a process (Harber and Weldon 1981, 1982). 

Singh and Pathria (1984) considered the BEC of a relativistic conformally coupled 
scalar field in the Einstein universe and found qualitative and quantitative agreement 
with Altaie (1978). Parker and Zhang (1991) considered the ultra-relativistic BEC 
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of the minimally coupled scalar field in an Einstein universe in the limit of high 
temperatures. They showed, among other things that ultra-relativistic BEC can occur 
at very high temperature and densities in the Einstein universe, and by implication 
in the early stages of a dynamically changing universe. Parker and Zhang (1993) also 
showed that the Bose -Einstein condensate could act as a source for inflation leading 
to a de Sitter type universe. 

Trucks (1998) repeated the calculations of Singh and Pathria but this time for 
the minimally coupled scalar field obtaining similar results. In fact performing the 
calculations for the minimally coupled scalar field amounts to substituting m = 
(m 2 + 1/a 2 ) 1 ^ 2 where m is the mass of the conformally coupled field. But as the 
calculations were effectively considered in the large radius region, the results come 
out to be identical with those of Singh and Pathria. 

The importance of the study of BEC in curved spaces stems from the interest 
in understanding the thermodynamics of the very early universe and that such a 
phenomenon may shed some light on the problem of mass generation in the very 
early universe. Indeed the investigations of Toms (1992, 1993) have shown that a 
kind of symmetry breaking is possible. 

In this paper we will consider the onset of BEC of the relativistic spin 1 particle- 
antiparticle system in an Einstein universe, a state which is surly relevant for the 
early stages of the universe at a point when the electromagnetic interactions decouple 
from the weak interactions. We will carry out the calculations in a similar fashion to 
that of Singh and Pathria and compare the results with our earlier non-relativistic 
case. Throughout this work we adopt the absolute system of units in which c = G = 
k = h = 1. 

2 The Charge Density 

We consider an ideal relativistic Bose gas of spin 1 confined to the background ge- 
ometry of the spatial section S* 3 of an Einstein universe with radius a. Since we are 
considering a relativistic system, it is necessary to take in consideration the possi- 
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bility of pair production. The system is taken to be formed of Ni particles and N2 
antiparticles. The total number of particles Q = N\ — N 2 is assumed to be conserved, 
though Ni and N 2 may change. A chemical potential u (is assigned for the particles 
and — fi for the antiparticles (see Harber and Weldon 1981). Accordingly the particle 
and antiparticle distributions are given by, 



iV 1 = £ci„[e^- 



-1 



N 2 = Y,d n [e^ 



+a0 _ 1 



where (3 = 1/T, e n are the eigen energies, and d n is the degeneracy of the nth level. 
The equation of motion of the spin 1 field in an Einstein universe was considered by 
Schrodinger (1938) and the solution yields the following energy spectrum, 



£n = - ( n2 + 
a v 



with degeneracy, 



d n = 2{n 2 - 1) , 
The charge density q is then found to be, 

Q 
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where V = 2n 2 a 3 is the volume of the spatial section of the Einstein universe, and 
(3' = l/Ta. In order to carry out the summation over n in (|j) we apply the Poisson 
summation formula (see Titchmarsh 1948), 
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These integrals can be easily performed using (Gradshteyn and Ryzhik 1980) 
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where K v are the modified Bessel functions of the second kind. 

Evaluating Iq and Ij and then substituting in (Q) we find that the charge density 
can be written as, 
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where 



z = yjp+f* and f = (27ra//3) j . 
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The first term in (fL0|) arises because the n = term in the summation over n in (^) 
is non-zero. The second term arises because of the spin-curvature coupling, which 
is defined mathematically by the coefficient a n of the Schwinger-De Witt expansion 
(see, De Witt 1965). The third term is just twice that of the scalar case considered by 
Singh and Pathria (1984). Both the first and the second terms disappear in the limit 
a — > oo. The bulk term is obtained in this limit assuming that q remains constant, 
which is the known thermodynamic limit. This reduces to the j = contribution of 
the last term in (|TU|) which gives, 
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This is just twice the value obtained for the scalar would be expected. The 

summation over j in (|10"D can be performed using the Poisson summation formula (|5|) 
again, where we obtain, 
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where // = /x + 27iip. The integrals in the above equation can be evaluated exactly 
using (see Singh and Pathria 1984), 
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Then using (see Gradshteyn and Ryzhik 1980), 
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where r = \Jm 2 — /j' 2 . 

The summation over Z can be easily done, so the exact expression for the charge 
density becomes, 
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Note that this form of the charge density is exact and no approximation whatsoever 
has been made through the calculation. 



3 Bose-Einstein Condensation 



We will adopt the microscopic criteria for marking the onset of the condensation 
(Altaie 1978), according to which the condensation region is defined such that a 
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large number of particles is found occupying the ground state. This implies that the 
chemical potential \l of the system approaches the minimum single particle energy 62, 
not the single rest mass energy considered by Singh and Pathria (1984). This can be 
observed directly from ([!]). However the consideration of the criteria that fi — > m is 
justified only for the minimally coupled scalar field case where the minimum energy 
is m (see Parker and Zhang 1991). In our case the chemical potential must satisfy 
the condition that 
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whereas in the conformally coupled spin case the condition is, 
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However, if m? is much larger than I /a 2 then it will be justified to take the limit 
fi — > m, but this approximation may impose certain restrictions on the range of 
the region under consideration. We will follow Singh and Pathria and adopt such 
approximation in this work. In such a case the main contribution of the summation 



over p in (|18"D comes from the p = term. Other terms are of order of e 
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0(e q / At ) where At = y2irf3/m is the mean thermal wavelength of the particle. The 
bulk term will reduce to (Singh and Pandita 1983), 
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Therefore we can write 
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If we define the thermogeometric parameter y as 
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then equation ( p2]) becomes 
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From this equation the behavior of the thermogeometric parameter y in the conden- 
sation region can be determined. It is clear that the second term in ( |24] ) defines the 
finite-size and spin-curvature effects. 

4 The Condensate Fraction 

The growth of the condensate fraction is studied here in comparison with the bulk 
case. This will show the finite-size and the spin-curvature effect in the range consid- 
ered, i.e., ma ^> 1. The charge density in the ground state is obtained if we substitute 
n — 2 in dT|). This gives 
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In the condensation region \i — > e 2 , so that the main contribution to the charge 
density in the ground state comes from the first term in the square bracket. This 
means that 

3 , x 



From (p3|) as /i — > m , we have 
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On the same footing we can expand e 2 as 
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So that (E6|) becomes 



6777 

-Mir*) ' (29) 
This means that the macroscopic growth of the condensate will occur only when 

y 2 — > — 4?t 2 , (i.e., y —>■ 27u). In order to see how this condensate compares with the 

bulk case we use the expansion: 
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Substituting this in (|4|) and using ( p9~D we get, 
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For the bulk system (a — > oo) there exists a critical temperature, T = T c given by 
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This condition can be written as 
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Thus we can write the bulk condensate density as 
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For the case of the finite system under consideration the condensate density is given 
by 
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where the y dependence on T in the condensation region can now be determined 
explicitly by 
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This equation can be solved numerically for the values of y at different temperatures. 
However, we notice that at T = T c (W(/3,m) = W(/3 c ,m)), the value of y can be 
obtained by solving the equation, 
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which has a solution at y c = 4.859i. 



5 Non-relativistic and Ultra-relativistic Limits 



The non-relativistic limit is obtained by setting (3m ^> 1. In this case we can use the 
asymptotic expansion of the Bessel functions of the second kind for large argument, 
where we have (see Abramowitz and Stegun 1968), 
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Therefore from (R6j) we have for the bulk system, 
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In the non-relativistic limit (|38|) reduces to 



2 ( - + - ] cothy - -J = x l/2 
Vtt y y 2 



1 — x 



-3/2 



j[2C(3/2)] 



2/3 



(40) 



(41) 



(42) 



(43) 



(44) 



where x = T/T c and I is the mean inter-particle distance. Note that equation ( fP| ) is 
the corrected version of equation (62) of our earlier paper (Altaie 1978). 
In the ultra-relativistic limit (3m 1, therefore, we use the expansion, 
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Accordingly the ultra-relativistic behavior of the bulk charge density given in (|36|) 
will be 
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However, substituting ( Bq ) in 
cles is 



the critical temperature for the bulk spin 1 parti- 

T -=JE- (48) 
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Figure 1: The quantity v = \fy 2 + 4n 2 , where y is the thermogeometric parameter, 
is plotted as a function of the scaled temperature x = T/T c for two different values 
of the radius a. The solid line is for the larger value of a. The vertical dotted line at 
x — 1 is marked for reference. 

This is the analogous result to that of the minimally coupled scalar field. This shows 
that the treatment of the problem using the large radius approximation as adopted by 
Singh and Pathria and in this work is equivalent to the high temperature expansion 
of Parker and Zhang (1991). This can be understood in the light of the fact that in 
the ultra-relativistic regime T oc a -1 , so that aT =constant. In the ultra-relativistic 
regime the thermogeometric parameter y behaves according to 

(i + i)c^ v = m( x -l) =m M( x -l). (49) 

\7r y J ml c \ xj V 3m V x) 

This equation can be solved for given values of a, q, and m. For illustrative purposes, 
figure 1 shows the behavior of the thermogeometric parameter y (drawn as v — 
\Jy 2 + 47r 2 versus the scale temperature x = T/T c ). It is clear that as a increases the 



10 



quantity y/y 2 + 47r 2 tends to a step function. 



6 The Critical Radius 



Although it is known that the universe as a whole is neutral or almost neutral, the 
observational limit on the net average charge density do not exclude the possibility 
that the charge density in the early universe was sufficient to produce relativistic 
Bose-Einstein condensation. Following the assumption that the charge is conserved 
we can write 



where qi and q v are the initial and the present charge densities, respectively. If we 
consider Ta = C and assume that C remains constant throughout the development 
of the universe then from (|48|) we can allocate a critical radius for the universe a c 
below which the gas will be always in the condensate state. Parker and Zhang (1991) 
have already noted this. However for the spin 1 case this critical radius will be given 



The estimated upper bound on the net average charge density of the universe at 
present is q < 10~ 24 cm -3 (see Dolgov and Zeldovich 1981). If this upper bound 
on the present charge density is adopted, then an upper bound on the radius of the 
universe at which condensate starts can be deuced. Using a value of 10 28 for C (Trucks 
1998) and a p = 10 28 cm, the critical radius for the onset of the condensation of the 
heavy gauge bosons W can be calculated. This gives a c < 10 -11 cm. 

7 Discussion 

In the previous sections we have investigated the behavior of an ideal relativistic spin-1 
gas confined to the background geometry of an Einstein universe. The extension of the 
problem into the Robertson- Walker spacetime is straight forward and results will be 
similar to the ones obtained here apart from the fact that in the Robertson- Walker 
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spacetime the radius a is time-dependent, so a — > a(t). On the other hand more 
accurate calculations may be needed in order to understand the behavior of bosons 
at very early stages of the development of the universe. In this context the problem 
of photon condensation will surely be of great interest, however some technical and 
conceptual problems hinders obtaining exact results. 
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